Binary search trees for generalized measurement 
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Generalized quantum measurements (POVMs or POMs) are important for optimally extracting 
information for quantum communication and computation. The standard realization via the Neu- 
mark extension requires extensive resources in the form of operations in an extended Hilbert space. 
For an arbitrary measurement, we show how to construct a binary search tree with a depth loga- 
rithmic in the number of possible outcomes. This could be implemented experimentally by coupling 
the measured quantum system to a probe qubit which is measured, and then iterating. 
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I. INTRODUCTION 



A crucial element of quantum information processing 
(QIP) and communication is measurement of a quantum 
system to access its information content, hence determin- 
ing optimal measurements is important. As QIP steps 
out from the pages of theory and into the laboratory, 
one has to find implementations given the usual limited 
resources of the real world. The most general measure- 
ment one can perform on a quantum system is given by 
a positive operator valued measure (POVM) [l| which 
can be considered as a projective measurement on an ex- 
tended Hilbert space of which the original state resides 
on a (proper) sub-space. Its realization via the Neumark 
extension 0, Q requires, broadly speaking, that the ex- 
tended Hilbert space should have as many dimensions 
as there are possible outcomes of the POVM. This has 
been described for atoms or ions [H, [6fl and for linear 
optics @, and POVMs have been realised on optically 
encoded quantum information d, [9], [T3, H3 ■ 

For many physical systems, however, it is difficult or 
impossible to find enough extra dimensions, let alone per- 
form operations on the extended system, hence a more 
efficient method is desirable. Recently, a method was 
discovered which allows an arbitrary POVM to be per- 
formed by adding only a single extra dimension to a sys- 
tem, essentially checking the measurement outcomes one 
by one However, when the number of possible mea- 
surement outcomes becomes large, more time efficient 
measurement strategies, also carrying a minimal dimen- 
sional overhead, would be useful. It is clear that a se- 
quence of partial conditional measurements implements 
a final effective POVM with many elements [1, Here 
we show, given any POVM, how to construct a suitable 
binary search tree of two-outcome POVMs by coupling 
the original system with a single qubit [22j . This way, a 
measurement with N = 2* outcomes can be implemented 
in t steps, resulting in a significant speedup. Existing ex- 
perimental realizations could easily be adapted to this 
method flq]. 



II. GENERALIZED MEASUREMENT 



A quantum measurement is often considered to be 
a projection in a complete basis of the d-dimensional 
Hilbert space. However, many experimental measure- 
ments are not well described by this. More generally, we 
only require of a measurement that the outcome prob- 
abilities are positive and sum to one, and satisfy con- 
vex linearity over mixtures of states. This leads to the 
framework of generalized quantum measurements, where 
a measurement is represented by a set of positive oper- 
ators {Mj}jLu (ip\Mj\i/j) > which sum to unity, 



Each outcome j is associated with an opera- 



tor Mj, and occurs with probability pj = Tr[Mjp], where 
p is the measured state. Hence a generalized measure- 
ment is usually called a positive operator valued measure 
(POVM) or probability operator measure (POM). 

The Neumark extension provides a way of performing 
any POVM via projective measurements, albeit in an ex- 
tended space. Without loss of generality, assume that 
each measurement operator Mj is proportional to a one- 
dimensional projector Mj = \il>j){t^j\ where is not 
neccessarily normalized [23] . If N is the number of out- 
comes and d the dimension of the Hilbert space, then 
N > d will hold. If we form the d x N rectangular array 
(M)jk — (k\ipj), where {\k)} is the computational basis, 
then the completeness relation implies that the columns 
of (M) are orthonormal V-dimensional vectors. Hence 
(M) can be completed to an V-dimensional unitary ma- 
trix Um whose j th row represents a state \ipj Xt ) in an 
iV-dimensional extended Hilbert space. The normalized 
projector \i/jj Xt ){i/jj Xt \ corresponds to outcome j for the 
original system. This procedure corresponds to applying 
U M to the extended Hilbert space in which the original 
system is embedded, and then making a projective mea- 
surement in the computational basis. If N is large, it 
may be infeasible to manipulate the required extended 
quantum system all at once, and we will therefore look 
at a way to reduce the number of ancillary dimensions 
by making sequential measurements. 
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III. SEQUENTIAL MEASUREMENT 

The {Mj} are sufficient to determine the measure- 
ment probabilities, but the post-measurement state is 
not uniquely defined. However, for any realization we 
can find Kraus operators {mj}, where Mj — m^rrij Vj, 
which tell us how the quantum state is affected 0, Qj} • 
If outcome j is obtained, then the quantum state p is 
transformed to pj = mjpm^/Tt\mjpm-]. A subsequent 
measurement, in general depending on the outcome j, 
acts on this transformed state. 

A series of measurements can be viewed as an ef- 
fective single generalized measurement, the sequence of 
outcomes determining the cumulative measurement op- 
erator. If the sequence ji,j2,~-,jt of outcomes have 
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, Mj t } and {to*^} as corresponding mea- 



surement and Kraus operators, then the final effective 
measurement operator and Kraus operator are given by 
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Here we assume that the measurement operators and the 
Kraus operators are d x d operators, i.e. the measure- 
ments map the system to a Hilbert space of the same 
dimension [24|. Hence, a sequence of measurements re- 
quires non-destructive measurement, e.g. indirect mea- 
surement of a system by measuring a probe after it has 
interacted with the system. 

A series of binary outcome measurements is shown 
in Figure [TJ The simplest probe is a two-level system 
(qubit), giving a binary measurement, a d-level probe al- 
lowing a d-outcome measurement. A unitary operator 
couples the probe with the system, e.g. via a coupling 
Hamiltonian over a set period. This in general entan- 
gles the state of the probe with the state of the system. 
Measuring the probe performs an indirect measurement 
of the system. From the Stinespring dilation [ill ], this 
effectively implements a completely positive map with 
Kraus operators given by bj = (j|?7|0), where {\j}} is 
the computational basis of the probe. Outcome j corre- 
sponds to the measurement operator Bj = b .bj and the 
conditional post-measurement state is 



Pi 



Tr[b jP b\] 



(2) 



By choosing suitable unitaries, any binary outcome 
POVM can be implemented at each stage. 

Conditioned on the result of the hrst measurement, a 
second measurement is performed, a third, and so on 
(Fig. [T]). This builds up a binary measurement tree 
with each pair of branches representing a different binary 
POVM, depending on the previous results. Each node 
represents the effective measurement operator (given by 
Eq. Ilbp obtained at that point. Hence, after t measure- 
ments, the effective POVM may have as many as N = d 1 
elements at the lowest level for a (i-level probe. 
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FIG. 1: Sequential POVM. A probe system in the state |0) in- 
teracts with the measured system p via unitaries. The probe 
is measured in the computational basis. Subsequent measure- 
ments, determined by U x , are conditional upon the preceding 
results i,j,k,.... The probe is reset before each measurement. 




FIG. 2: Four-outcome POVM. The effective measurement op- 
erators {Mij} are given by the B measurements at each step. 
We want to determine what binary measurements B are re- 
quired to implement {Mij}. 



IV. BINARY MEASUREMENT TREES 

It is easy to build up POVMs with many elements from 
a binary measurement tree. However, given an arbitrary 
POVM with elements Mj, constructing such a measure- 
ment tree which implements it is not so obvious. Here 
we show how it may be done. 

It is instructive to look at the simplest non-trivial bi- 
nary POVM tree with t = 2 (Fig. E}. Let B { and By- 
denote the binary measurement operators performed at 
the first and second stage, and Mj,Mjj denote the cu- 
mulative measurement operators. The following should 
hold, where j, k = 0, 1: 



Mj = 
I = 



Bo - 



Mj = M j0 



m-ik 



Bi — Bj 



B 



biking 



(3a) 
(3b) 
(3c) 
(3d) 



Let us take bo = too, b\ = toi and use the ansatz 6y = 
mijb~ l where the Moore-Penrose pseudo-inverse A^ 1 of 



an operator A is uniquely defined by [191 ] 



AA~ 1 A 
A^AA- 1 
AA- 1 
A- X A 



A, 

(AA-y, 



(4a) 
(4b) 
(4c) 
(4d) 



We shall prove that the bij so constructed correspond to 
POVM operators and solve the task. 

First, since Mj is a positive operator, Mj = 
Sfc=i ^fe| e ife)( e *|j with positive eigenvalues \k and cor- 
responding eigenvectors \en-); r is the rank of Mj. The 
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FIG. 3: Binary POVM Tree. An arbitrary POVM is given by 
operators Mj — rmmj. A sequence of binary outcomes with 
measurement operators J3 M , B v ,...,£?{ leads to the measure- 
ment Mj = m^rrij where mj = 6j . . . b^b^. Half of the possible 
results in the branches below are eliminated at each step. 



Mij are positive operators and JV My = Mi, so the null 
space of Mi is contained in the null spaces of My , hence 
Mij = Y^ r k,i=i4>ij,ki\ e ik){eii\ for some </>y,fcz- Similarly, 
m ij = J2li=i <Pij,ki\ e ik)( e ii\ for some <p ijtk i- 

We can expand b t = m l = Y] fe \Z%ikVi\e ik ) {e ik \ for 
some unitary V l , similarly 6" 1 = J2k l/\/^l e ;fc)(eifc|V/. 
Hence, we can see that 




{=1 



In general, completeness of {i?y}j requires us to modify 
our original ansatz by adding an extra operator, 

bij = rriijb' 1 + a 2 g { , (5) 



ond level, from the previous section we have 

&12...W/4 = m 12...N/J> 1 2...N/2 + 912...N/2 

&JV/4+1..JV/2 — m A r /4+l...A , /2^i 2 1 ..Ar/2 + 512...JV/2 

&i\r/2+1...3jV/4 = TO Af/2+l. .. 3iV/4^jv/2+l. .. JV + 9N/2+1. ..N 

&3JV/4+X...JV = m 3A r /4+l...Arfc A r/ 2 +l...JV + 9N/2+\. ..N 

where we have absorbed the normalization of the oper- 
ators. At subsequent levels, we can express the required 
binary POVMs as 

b Xa = m Xa m~l Xb +g XaXi (6a) 
b Xb = m Xb m-l Xb + g XaXb (6b) 

where x a Xb is the concatenation of the strings x a and Xb- 
At the last level bi = mpjj 1 + <?i2 and 62 = JT^JTiTjg +<7i2- 
Note that the unitary freedom m x — > Vjirij leaves the 
observed probabilities invariant but simply rotates the 
post-selected states after each measurement. 

For an N element POVM, we need only a probe qubit 
and [log 2 JV] rounds of binary measurements. Let us 
determine the number of operations required to imple- 
ment this measurement compared to other methods. For 
a measurement with N outcomes on a d-dimensional 
quantum system, the standard Neumark extension re- 
quires a N x JV unitary transform. This can be real- 
ized with N(N ~ l)/2 operations between pairs of basis 
states [IS], followed by a projective measurement in the 
JV-dimensional space. The realization using just a single 
extra degree of freedom [1] requires a (d+1) X (d+ 1) uni- 
tary transform to be implemented a maximum of N — d 
times, giving in total a maximum of (N — d)(d + l)d/2 
operations [25[. The binary search requires a 2d x 2d 
transform to be implemented [log 2 JV] times, that is, 
[log 2 JV] d(2d — 1) pairwise interactions, a significant 
speedup if JV is large. 



where gi — Xy= r +i l e y)( e y'l^ ^ s an isometry on the 
null space of b\ and the coefficients satisfy JV \ a j\ 2 = 1- 
We have defined gi so that gibi = 0. With this slight 
modification, it is easy to show that J^ fc Bjk = I. 

For a general POVM {M k } with K elements, we first 
pad the set with null operators until it contains JV el- 
ements for JV = 2*,t = |~log 2 K~\ (Fig. [3]). In a conve- 
nient change of notation, the cumulative POVM at the 
j th level consists of 2 J operators M x where x is a se- 
quence of 2*~ J numbers indicating which of the possi- 

ble outcomes M x — JJ i=1 M Xi sit in the corresponding 
branches below. A binary POVM {B Xa , B Xb ] splits each 
node into two possible branches, each containing half of 
the remaining outcomes. We now determine the binary 
POVMs B which take us from a higher to lower branch. 

At the first level, B X2 ...n/z = 5Zi=i M i = M2...W/2 
and B N/2 ...N = Z)i= n /2+i M * = ^JV/2+i...jv- At the sec- 



V. EXAMPLE: TETRAD MEASUREMENT 

As an example of the method, consider the symmetric 
informationally complete POVM of a single qubit, the so- 
called tetrad measurement, with measurement operators 

Ma = 



l^iX^jl S iven by [21 1 
IV>o> 



>' 



IV'i) 
IV*) 



-LdO+^e^ll)), 

-L(|o) + ^ e ^|i)), 

^(|0) + v^|l)). 



(7a) 
(7b) 
(7c) 
(7d) 



Although the tetrad POVM can be performed in one pro- 
jective step with the addition of just one extra qubit, it 
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FIG. 4: Tetrad POVM. a) The first binary measurement is 
a partial filtering with operators {Mo3,Afi2}. b) The sec- 
ond (projective) measurements depend on the outcome of the 
first measurement. If M03 (M12) was obtained, then {-Bo, B3} 
({Bo , B3}) is measured. The second binary measurements are 
projective measurements in the plane perpedicular to the di- 
rections of the first measurement and the Bj lie in the direc- 
tion of the projection of the Mj upon this plane. 



is instructive to demonstrate the binary tree approach 
using this example. 

At the first stage, we are free to choose which two final 
measurement operators to group together, for instance, 



M 03 = B03 = Mo + M 3 = 



M12 = B 12 = Mi + M 2 = - 




(8a) 
(8b) 



We are also free to choose the Kraus operators m x = 
y/M x using the singular value decomposition. For exam- 
ple, 



m 03 = v A +l e +)( e +l + v A -l e -)( e -l, ( 9a ) 

m u = v /A _ |e+)(e + | + v/A7|e_)(e_|, (9b) 



where the eigenvalues and eigenvectors are 
A± = (l±V3)/2, 
\e±) = 



(10a) 

(±V3±V5|0) + VW\/3|1))/V6. (10b) 



Although Mq3 and M12 share their eigenbases, we 
need a full 4x4 Neumark extension binary POVM so 
that the post-measurement state is ready for the next 
stage. We couple the system via U to a auxiliary probe 
qubit prepared in the state |0) a . Then, projecting the 
probe onto states |0) a and |l) a corresponds to operations 
Q (0|i7|0) Q = m 03 and a (l|t/|0) a = m i2 on the system. A 
suitable coupling U is constructed by making a 4 x 4 
Neumark extension of the two-column matrix with its 
first two rows given by mo3, and last two rows by TO12. 
In the basis {|e±)|j) a }, one possible U is 



U 



A+ 

o_ yxr 

A_ 





A_ 











(11) 



In this example, the positive operators rrijk are invert- 
ible so the bj for the next step are easily obtained as 



b 



V^\/Mo3 , h 
b 2 = V^VM^ ~ , 63 



r M[y/Mv2 (12a) 
/MgVMoT 1 , (12b) 



which gives 

Bo = 

B x = 
B 2 = 

B, = 



1 / 1 



-1 



-1 



1 ( 1 -i 

i 1 



b\b 2 



1 ( 1 



-i 1 



1 / 1 



1 



1 

V^5 



1 - 



(13) 

(14) 
(15) 

(16) 



The bj are rank one operators but are not Hermitian. We 
can visualize the sequence of measurements on the Bloch 
ball (Fig. ID). 



VI. CONCLUSION 

In conclusion, we provide a constructive proof of the 
universality of sequential two-outcome POVMs. We show 
how to construct binary measurement trees to imple- 
ment any generalized measurement through a sequence 
of indirect binary POVMs requiring only an extra aux- 
iliary qubit. This avoids having to manipulate extended 
Hilbert spaces (larger than twice the dimension of the 
measured system) and reduces the number of required 
operations when the number of outcomes becomes large. 
The number of steps is logarithmic in the number of mea- 
surement outcomes. The required interaction to perform 
binary POVMs exists in physical systems such as cavity 
quantum electrodynamics (CQED) [l6| where the state 
of a field can be probed by an atom-cavity interaction 
and the atom measured. So far, projective measurements 
have been performed with a fixed interaction and mea- 
surement, but it should be possible with feed- forward 
and suitable control fields to implement a full POVM 
measurement as described here. 
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